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We report the observation of a new Daj meson produced in B+ D°Dsj -» D° D° K+ . This 
state has a mass of M = 2708 ± 9liJ MeV/c^ a width T = 108 ± 2Ztfi MeV/c^ and a 1~ spin- 
parity. The statistical significance of this observation is 8.4(j. The results are based on an analysis 
of 449 million BB events collected at the T(4S) resonance with the Belle detector at the KEKB 
asymmetric-energy e^e^ collider. 

PACS numbers: 14.40.Lb, 13.25.Hw, 13.20.Fc 



At the level of quark diagrams, the decay B DDK 
proceeds dominantly via the CKM-favored b cW^ 
CCS transition. The transition amplitudes can be cate- 
gorized as due to either external W- or internal (color- 
suppressed) M^-emission diagrams. The decay 
D'>D''K+ m can proceed through both types of dia- 
grams; thus it is promising for searches for new cs states 
as well as for some cc states lying above D^D^ thresh- 
old. The unexpected discoveries of the -D*q(2317) and 
1^51(2460) mesons show that our understanding of cs 
spectroscopy might be incomplete, while experimental 
data on cc states with decay channels open to D^*^D'^*^ 
are scarce. 

The decays B — > DDK have been previously studied 
with a small data sample at LEP 0] and more recently 
a larger statistics exploratory study was performed by 
BaBar In this letter we report the first study of the 
Dalitz plot of B+^D°D°K+. 

The study is performed using data collected with the 
Belle detector at the KEKB asymmetric-energy e^e~ 
(3.5 on 8 GeV) collider operating at the T(4S') reso- 
nance {y/s ~ 10.58 GeV). The data sample corresponds 
to the integrated luminosity of 414 fb~^ and contains 
449 million BB pairs. The Belle detector is a large-solid- 
angle magnetic spectrometer that is described in detail 
elsewhere 

Well measured charged tracks are identified by combin- 
ing information from time-of- flight, Cherenkov and ion- 
isation detectors. Requirements on the particle identifi- 
cation variable are imposed that identify a charged kaon 
with 90% efficiency, a charged pion with almost 100% ef- 
ficiency and have less than 10% K'^tt misidentification 



probability. Any track that is positively identified as an 
electron is rejected. 

Candidate Kg ^ tt'^tt' decays arc identified by a dis- 
placed secondary vertex, a two-pion momentum vector 
that is consistent with a Kg originating from the IP and 
a TT+TT^ invariant mass within ±15 MeV/c^ (i3cr) of 
the nominal Kg mass. Candidate tt*^ mesons are recon- 
structed from pairs of identified photons, each with a 
minimum energy of 50 MeV, that have an invariant mass 
within ±15 MeV/c^ (±3cr) of the 7r° mass. 

D*^ mesons are reconstructed in the K^tt'^, 
A'~7r+7r+7r~, A'~7r+7r°, AT^tt+tt" and K^K+ de- 
cay modes. We preselect D° candidates using a signal 
window ±30 MeV/ (±5cr) around the nominal D° me- 
son mass for all decay modes except for D° — > A'~7r+7r°, 
where a ±50 McV/c^ (iScr) signal window is used. 
Mass- and vertex-constrained fits are applied to D^ 
candidates to improve their momentum resolution. 

To suppress the continuum background (e+e^ qq, 
q = u,d,s,c) we require the ratio of the second to the 
zeroth Fox- Wolfram moments Q to be less than 0.3. 

We form Z)°I?°A'+ combinations using D^ candidates 
with momenta in the T(4S') rest frame (cms) kinemat- 
ically allowed in B+ D^D'^K^. The momenta of 
the secondaries from a B meson candidate decay are 
refitted to a common vertex with an interaction point 
(IP) constraint that takes into account the B meson de- 
cay length. The B meson candidates are identified by 
their cms energy difference, A A = S^A; — Abeam, and 
their beam constrained mass, Mbc — V-E-beam ^ (^iPO^, 
where Abeam = -v/s/S is the beam energy in the cms and 
Pi and Ei are the three-momenta and energies of the B 
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candidate's decay products. We select B candidates with 
Mbc > 5.2 GeV/c^ and -0.4 GeV < A£; < 0.3 GcV. 
Exclusively reconstructed i?+ D^D^K^ signal events 
have an distribution that peaks at the nominal B 
meson mass; the AE distribution peaks at zero. 

We employ a discriminator (likelihood ratio) based on 
the meson signal significance to select the unique 
B candidate in the event, defined as: C'lZ{M£io) = 

S(M ^o)+B{M ) ' '^'^^^^ ^ ^^'^ ^ ^""^ signal and the 
background likelihoods that depend on the D° candi- 
date's invariant mass (-M^jo). This discriminator is de- 
termined from fits to the AIjjo distributions for each 1)° 
decay mode separately, using a data sample enriched 
in i?+ D^D^K^ decays. In these fits S and B arc 
parametrized respectively by a double gaussian and a lin- 
ear functions. For events with multiple B^ D^D^K^ 
candidates, the product CRb = CR[MQa) x CR[Mjja) 
is calculated and the candidate with the largest CR-b is 
accepted. The CTZb discriminator is also used to sup- 
press combinatorial backgrounds to B'^ — *■ D^D^K'^ and 
to enhance the signal purity. Monte Carlo (MC) studies 
showed that this selection, which does not rely on the AE 
and Mbc values used in the B signal definition, does not 
introduce biases. The fraction of events with multiple 
candidates is 45% and the average candidate multiplicity 
is 2.2. The main background is due to BB production, 
with multiple candidates originating from wrong pairing 
of Z?°'s or swapped kaons. 

The AE and Mbc distributions for the B+-^D°D°K+ 
decay candidates, selected with CTZb > 0.01 require- 
ment, are shown in Fig. [l] where the AE distribution 
is shown for | Mbc - tob |< 3a m^,^ (""Mbo =2.7 MeV/c^, 
tob is the nominal B meson mass) and the Afbc distri- 
bution is shown for | AE |< 3a ae (cab — 6.6 MeV). 
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FIG. 1: AE vs. Mbc (a), AE (b) and Afbc (c) distributions 
for B+ ^D°D°K+. Black (wliite) ellipses in (a) enclose the 
signal (sideband) regions described in the text. 

From a study of the Mbc and AE background distribu- 
tions in large MC samples of generic BB and qq events as 
well as ZJ^-mass sidebands in data, we find no significant 
peaking background. 

To extract the signal yield, we perform two- 
dimensional (2D) extended unbinned maximum- 
likelihood fits to AE and Mbc- The probability density 



functions (PDFs) for the Mbc and AE signals are 
Gaussians. The background PDF for Mbc is represented 
by a phenomenological function 0] with a phase-space- 
like behaviour near the kinematic boundary; the AE 
background is parameterized by a second-order polyno- 
mial. The likelihood function is maximized with free 
parameters for the signal yield, the Gaussian means and 
widths, and four parameters that describe shapes of the 
background distributions. From the fit, we obtain a 
signal yield of iVsig = 399 ± 40 events. The results of the 
fit arc superimposed on the AE and Mbc projections 
shown in Fig. [Hb)-(c). 

We determine the branching fraction from the relation: 
B{B+ D°D°K+) = - — ^ ^l^n :.„,r.n where 



Cij are efficiencies for the decay channels i and j. Nbb 
is the number of analyzed BB pairs, Nbb = 449 x 10^, 
and Nb+b- = ^50^0 is assumed. The efficiencies are 
determined by MC using a model that reproduces the 
observed Dalitz plot features (discussed below). The sum 
in the denominator of the above relation is 4.0 x 10~*. 
Wc obtain B{B+~*D°D°K+) = {22.2±2.2tli) x 10"", 
where the first error is statistical and the second is sys- 
tematic. The latter includes contributions due to un- 
certainties in the efficiency determination (tracking and 
particle identification efficiency, data-MC differences in 
Ai?, Afbc signal shapes), the CTZb selection, the back- 
ground parameterization, the MC model used in the ef- 
ficiency calculation, the intermediate D'^ branching frac- 
tions and Nptpi. This result supersedes our previous de- 
termination 111, which assumed a phase space model in 



the efficiency determination. 

The main features of the data can be seen in the Dalitz 
plot M'^{D°D°) vs. M2(D°iC+) for events from a sig- 
nal region defined by the ellipse B? = (AE/aAE)'^ + 
{{Ahc~mB)/aM^^^, i?2< (1.5)2 shown in Fig. [^a). The 
three two-body invariant mass distributions are shown 
in Figs. [2I^b)-(d). The hatched histograms represent the 
background distributions obtained for events from an el- 
liptical strip surrounding the AE — Mbc signal region, 
defined by 6^ < i?^ < 10^. The background distributions 
are normalized to the number of background events un- 
der the signal peak (±1.5(j) as determined from the 2D 
AE and Mbc fit. The data are not efficiency corrected. 
The efficiency as a function of invariant mass is shown in 
Figs. [2l^b)-(d) as a continuous curve. 

A pronounced feature of the Dalitz plot is the accumu- 
lation of events in the region 16 GeV^/c-* < AP{D°D°) < 
18 GeVVc" and 7 GeY^ / < M^{D°K+) < 8 GeVVc'' 
possibly due to the overlap of a horizontal band that 
could be due to the -0(4160), ■0(4040) and a vertical band 
that cannot be attributed to any known cs state. A hor- 
izontal band at A'P{D°D°) ~ 14.2 GeVVc'* corresponds 
to 0(3770) production. 

The distributions in Fig. 2 are meant only to illustrate 
the features of the data. In the subsequent analyses a 
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FIG. 2: Dalitz plot (a) and projections for B+-*D°D°K+ in the 1.5a AE-My,^ signal region: M{D°D°) (b), M{D''K+) (c), 
M{D°K+) (d). Hatched histo grams represent background, red/solid curves show the efficiency. 



more robust procedure was used to obtain background- 
subtracted mass distributions. Wc determine the AE 
vs. Mbc distributions for events from mass bins of the 
Dalitz plot projection and fit the signal and background 
shapes to obtain B meson signal yield vs. invariant mass. 
20 MeV/c^ mass bins are used for the '0(3770), while 



The background-subtracted A/(Z3"Z)") spectrum 
(Fig. EJb)), for events satisfying cos^hci > 0, where 
^hci is the helicity angle between the D'^ momentum 
vector and the direction opposite the A'^ in the 
jjOjjQ j.gg^ frame, is used to estimate the possible 
■!/i(4160), V(4040) contribution to the enhancement 



50 MeV/c^ bins are used for the other studied samples. at M{D°K+) ~ 2.7 GeV/c^. The peak at thresh- 



The parameterizations and AE—Mtc ranges considered 
in these fits are the same as those used for the total B 
yield extraction. The widths and means of the Gaussians 
describing the signal are fixed at the values obtained for 
the total signal sample, while the signal yield and the 
background PDF's parameters are free parameters. 

The background-subtracted M{D"D") in the V(3770) 
signal region is shown in Fig. [3]Ja). The peak is fitted for 
M{D°D°) < 4 GeV/c2 with a Breit-Wigner (BW) plus a 
threshold function to describe a phase space component. 
The V'(3770) signal yield is 68 ± 15 events with a peak 
mass of 3776±5 MeV/c^, and a width of 27±10 MeV/c^, 
in agreement with the PDG averages. 
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FIG. 3: B meson signal yield vs.: (a) M(D°D°) in the 
V)(3770) region, (b) M{D°D°) for cos6lhoi > (c) M{D"K+) 
for M{D°D°) > 3.85 GeV/c^ Solid curves denote the 
fit results described in the text. The red/dotted curve in 
(a) shows the phase-space component, whereas in (c) the 
red/dotted curve is the sum of the three components: i/;(4160) 
reflection, phase-space and threshold components. 



old corresponds to ■0(3770), while the structure at 
4.0 4.2 GeV/c^ is conservatively assummed to be 
predominantly (0, ^) due to the V(4160). The distri- 
bution for M{D°D°) > 3.8 GeV/c^ is fitted with a BW 
with mass and width fixed at the nominal '!/'(4160) values 
(M = 4160, r = 80 MeV/c^ Q), yielding 24± 11 signal 
events. We use these ■0(4160) parameters to estimate the 
number of ^/;(4160) events in the backward helicity-angle 
hemisphere, in the region M{D'^K+) < 2.9 GeV/c^. 
Taking into account the efficiency we obtain a total of 
43 ± 20 ?/'(4160) events. 

Figure [3jc) shows the background-subtracted 
M{D°K+) distribution for events with M{D°D°) > 
3.85 GcV/c^. This requirement removes the -0(3770) 
reflection at high M{D°K+). The predicted 0(4160) 
reflection agrees well with the data in the high mass 
AI{D'^ K^) region but does not explain the large peak 
at M{D°K+) ~ 2.7 GeV/c^. We parameterize the 
observed excess of events with a BW and flt the 
M{D^K^) spectrum (Fig. [31(c)) with the ansatz of a 
new resonance, the -0(4160) reflection and a phase-space 
component with shapes determined by MC simulations. 
The efficiency variation is taken into account in the fit; 
the free parameters are the resonance yield, mass and 
width, and the phase-space component normalization. 
The fit has an acceptable overall but is unable 
to reproduce the events near the low-mass threshold 
seen in Fig. [3Kc). We used several phenomenological 
parameterizations (polynomials, a BW, an exponential) 
of the threshold enhancement in the fit to determine 
its influence on the BW parameters of the 2.7 GeV/ 
peak. The exponential form a x exp [— aAf ^(D°Ar+)] 
gives a good description of the mass spectrum, while 
adding only two free parameters. 

For the new resonance, which we henceforth denote 
as the D,5j(2700)"''; we obtain a signal yield of 182 ± 30 
events, a mass of M = 2708 ± 9 MeV/c^ and a width of 
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TABLE I: Resonance parameters and product branching frac- 
tions: B{B+ D°D,ji2700)+) x B(D,j(2700)+ ^ D°K+) 
and S(B+^i/»(3770)iC+) x B(i/>(3770)^D°D°). 



R 


D.J (2700)+ 


V'(3770) 


Nsig (Significance) 
M [MeV/c^] 
r [MeV/c^] 


182 ±30 (8.4cr) 
2708 ± 9t\l 
108 ± 23^3^ 


68 ±15 (5.5ct) 
3776 ± 5 ± 4 
27± 10 ± 5 


Product B [10"''] 


11.3 ± 2.2^2,1 


2.2 ±0.5 ±0.3 



r = 108±23 McV/c^. The threshold and the phase-space 
components from the fit are 58 ± 38 and 47 ± 26 events, 
respectively. The fit results are shown in Figs. |lja)-(c) 
as histograms overlaid on the measured mass spectra. 

The resonance parameters and product branching frac- 
tions are summarized in Table [J (the first error is statis- 
tical, the second is systematic). The systematic errors 
on the product branching fractions and the resonance 
parameters include contributions from uncertainties in 
the yields of the '0(4160) reflection (including the recent 
-0(4160) parameter determination [9|), the threshold pa- 
rameterization, sensitivities of parameters to the fit range 
and parameterization, uncertainties in the CTZb selec- 
tion, as well as uncertainties due to interference effects 
that were neglected. The systematics due to the lat- 
ter are determined from MC simulations of Dalitz plot 
densities with and without interference of contributing 
amplitudes, with each contributing resonance parame- 
terized by a BW form. The resonance parameters from 
Table U and the threshold enhancement parameters are 
used to determine the amplitudes. The effects of inter- 
ference of the -0(3770) with other states are found to be 
small and are neglected in the simulations. These MC 
samples, with maximal constructive and destructive in- 
terferences, were analysed ignoring interference effects. 
The differences between the obtained resonance parame- 
ters and the input values are taken as systematic errors. 

We study background-subtracted 0(3770) and 
Z)sj(2700)+ helicity angle distributions by selecting the 
respective invariant mass in the resonance region and ob- 
taining B meson signal yields in bins of cos ^hci from the 
2D fits to AE and Mbc. Here cos6'hei for ^(3770) is de- 
fined as before, whereas for Dsj{2700)^ it is the angle 
between the momentum vector and the direction op- 
posite the D° in the D^K^ rest frame. The obtained 
angular distributions are then corrected using bin-by-bin 
efficiencies. The expected refiections from 0(4160) and 
from the threshold component are subtracted from the 
£'sj(2700)+ angular distribution. Spin hypotheses for 
the resonances are tested by comparing predictions for 
the different hypotheses to the corrected angular distri- 
butions. The -0(3770) distribution (not shown) is well 
described by the J = 1 hypothesis {x^/ndf = 3.6/5). 
The i:»^j(2700)+ distribution (Fig. H^d)) favours J = 1 
(11/5); the J = (112/5) and J = 2 (146/5) assignments 
can be rejected. The J = 1 assignment and the observed 




2.6 2.8 _ 3 3.2 3.4 

M(D°K*) (GeV/c^) 



FIG. 4: B meson signal yield vs. : M{D'^D°) (a), 
M{D°K+) (b), and M{D°K+) (c). Histograms denote 
the contributions from: D.j(2700)+ (blue/grid), 0(3770) 
(green/horizontally striped), 0(4160) (yellow/light grey), 
threshold (red/vertically striped) and phase-space compo- 
nents (brown/dark grey). Histograms are superimposed ad- 
ditively. (d) Efficiency corrected Dsj(2700)+ helicity angle 
distribution. Curves show predictions for various spin hy- 
potheses: J — Q (green/dotted line), J=l (red/solid), J = 2 
(blue/dashed), where Legendre polynomial values are aver- 
aged over the bin width. 

decay to two pseudoscalar mesons imply parity P ~ —I. 

In summary, from a study of the Dalitz plot we find 
that the decay D'^D'^K^ proceeds dominantly 

via quasi- two-body channels: i3+ ^ Z)"Ds,/(2700)+ and 
B+ -0(377O)i^+. The observed rate for 0>(377O) pro- 
duction in B meson decays confirms our previous ob- 
servation [ll[. The Dsj(2700)"'" is a previously unob- 
served resonance in the D'^K^ system with a mass M = 
2708 ± 9liJ MeV/c2, width F = 108 ± 23^i^^ MeV/c^ 
and = 1~. The statistical significance of this ob- 
servation is 8.4(7. Based on its observed decay channel, 
we interpret the £'sj(2700)''' resonance as a cs meson. 
Potential model calculations jl^l predict a cs radially ex- 
cited 2^Si state with a mass 2710-2720 MeV/c^. From 
chiral symmetry considerations [l^l a 1+-1~ doublet of 
states has been predicted. If the 1+ state is identified 
as the Dsi(2536), the mass predicted for the 1^ state is 
M = 2721 ± 10 MeV/c^. Additional measurements of 
the meson properties are needed to distinguish between 
these two interpretations. 

It is not clear whether the structure at 2688 MeV/c^ 
observed recently 1^ in the DK system produced in con- 
tinuum could be due to the Dsj{21Q{))^ . The recently re- 
ported Dsj (2860) state [l3| is not seen in our data. This 
could indicate a high spin for this meson that suppresses 
its production in B decays. 
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